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Introduction 

A  knowledge  Is  often  required  of  the  highest  zenith  angle  for  which  a 
charged  particle  of  given  rigidity  has  access  to  an  earth  satellite  lying  in 
a  particular  orbital  configuration. 

Because  of  the  complexities  of  the  trajectories  by  which  below-llne- 
of  sight  access  to  a  satellite  occurs  it  has  proven  difficult  to  determine 
systematically  whether  absolute  limits  exist  beyond  which  access  Is  forbidden. 
The  general  approach  to  this  problem  has  been  to  directly  test  for  access  by 
carrying  out  a  large  volume  of  trajectory  traces.  This  method,  however.  Is 
not  completely  satisfactory  because  the  absence  of  detected  access  beyond 
a  given  given  angle  does  not  constitute  a  satisfactory  proof  that  access 
could  not  occur  (either  via  undetected  allowed  tractorles  that  exist  under 
the  geomagnetic  field  conditions  simulated,  or  under  even  slight  perturbations 
away  from  these  conditions). 

A  technique  has  been  developed,  based  on  the  trajectory  parameterization 
method  (Cooke,  1982),  which  allows  a  very  much  more  definitive  insight  into 
the  problem  of  establishing  the  high  zenith  angle  limits. 


Discussion 


The  new  Method  involves  the  examination  of  the  form  of  trajectories 
along  which  particles  approach  a  satellite,  and  consideration  of  whether 
these  orbits  can  or  cannot  provide  access  for  charged  particles  entering  from 
outside  the  geomagnetic  field.  (In  line  with  ordinary  trajectory  tracing 
practice  the  traces  are  carried  out  in  the  reverse,  negative-time,  direction 
-  out  away  from  the  satellite.) 

Before  discussing  the  technique  further  it  is  appropriate  to  briefly  re¬ 
view  the  question  of  access  to  a  point  in  an  axially  symmetric  representation 
of  the  geomagnetic  field  (as  Investigated  by  Stormer  (1930,  1955),  and  by 
Lemaitre  and  Vallarta  (1936a,  1936b).  In  a  simple,  axially  symmetric  field, 
three  types  of  solid  angle  regions  may  be  distinguished  relative  to  any 
point  in  the  field: 

a)  The  “main  cone**  -  the  region  into  which  access  is  directly  allowed 
via  simple  aperiodic  trajectories, 

b)  The  "Stormer  cone"  -  into  which  access  is  forbidden  to  particles 
from  outside  the  field,  In  this  region  trajectories  passing  through 
the  given  point  have  the  form  of  bound  periodic  orbits. 

c)  The  "penumbra"  -  the  solid  angle  region  lying  between  the  main  and 
Stormer  cones,  into  which  access  is  neither  expressly  allowed  nor 
forbidden  by  the  theories  of  Stormer,  and  Lemaitre  and  Vallarta.  The 
region  typically  consists  of  admixed  allowed  and  forbidden  subregions. 

In  the  real  field  (non  axially  symmetric)  the  situation  is  complicated 
by  the  fact  that  the  Stormer  cone  is  very  difficult  to  define,  because  the 
associated  trajectory  forms  do  not  correspond  to  perfect  bound  periodic 
orbits.  It  has  been  difficult  to  determine  the  properties  of  the  quasi 
bound  periodic  orbits,  and  to  decide  whether  access  from  outside  the  field 
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could  occur  along  such  trajectories.  This  problea  has  in  part  been  addressed 
In  the  present  study,  and  satisfactory  working  principles  have  been  developed 
for  dealing  with  the  question  of  possible  access. 

For  the  purpose  of  this  investigation  it  has  proved  practical  and  useful 
to  distinguish  four  different  kinds  of  solid  angle  regions  relative  to  a 
given  point  in  the  “real1*  field.  These  are  the  regions  within  which: 

(i)  access  is  forbidden  because  of  short  range  intersection  with  the 
atmosphere  (as  detected  in  the  outwards  trajectory  trace).  These 
trajectories  correspond  to  the  short  range  “shadow"  orbits  of 
Lamaitre  and  Vallarta. 

(ii)  access  from  outside  the  geomagnetic  field  either  does  or  could 
possibly  occur  along  aperiodic  or  unstable  quasl-periodlc  orbits, 

(iii)  entry  would  only  be  by  aperiodic  or  quasl-periodlc  orbits  that 
intersect  the  surface  of  the  earth. 

(iv)  approach  could  occur  via  quasl-periodlc  orbits  that  have  no  de¬ 
tected  intersection  with  the  atmosphere.  Such  trajectories  can 
extend  for  very  large  distances  in  the  field  without  any  detected 
field  escape  or  earth  intersection.  It  is  impractical  (or  even 
sometimes  impossible)  to  trace  these  to  some  definite  end,  so 
the  traces  are  terminated  at  a  nominated  path  length.  Thus  the 
matter  of  access  via  such  trajectories  remains  unresolved  (the 
question  of  the  precision  of  very  long  trajectory  traces  neces¬ 
sarily  enters  here,  but  will  not  be  discussed). 

Figures  1  through  4  illustrate  trajectories  associated  with  the  four 
kinds  of  region.  In  the  trajectory  parameterization  technique  used  in  this 
study  (in  which  trajectory  features  are  "recognized”  by  computer,  and  their 
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positions  Monitored  for  use  in,  for  exaaple,  napping  cut-off  structures, 
and  establishing  penuabral  structure  detail),  it  is  laportant  to  continue 
the  trajectory  traces  on  beyond  the  first  earth  intersection.  This  is  an 

artifact  necessary  to  the  operation  of  the  technique  (and  in  the  identifies-  ' 

tion  of  trajectory  fora),  and  is  not  held  to  represent  any  kind  of  reality. 

Thus,  in  Plgures  1  and  3,  the  trajectories  are  seen  to  pass  through  the 
surface  of  the  solid  earth  (where  the  location  of  the  surface  is  represented 
by  the  dashed  line). 

Prellalnary  investigations  have  been  carried  out,  using  the  trajectory 
paraaeterlzatlon  technique,  to  elucidate  the  fora  and  properties  of  the 
regions  containing  the  trajectories  of  the  four  kinds.  As  a  first  step  in 
this  study,  the  regions  containing  trajectories  of  type  (lv)  are  mapped, 
either  in  zenith-azimuth,  or  in  latitude-longitude,  space.  An  example  of 
zenith-azimuth  mapping  is  given  by  Cooke  (1982),  and  reproduced  in  part  as 
Fig.  5.  The  high  zenith  angle  Halts  are  set  by  the  structure  visible  in 

the  west  at  ~145°  zenith  angle.  This  boundary  is  formed  by  short  range 

shadow  orbits. 

Uhen  aapped  in  latitude  and  longitude  the  regions  typically  have  the 
form  shown  In  Figure  6.  At  the  given  zenith  and  azimuth  (150s  and  270° 
respectively)  access  by  positively  charged  7  GeV  particles  to  a  satellite  in 
a  400  km  geocentric  orbit  is  explicitly  forbidden  at  all  points  on  the 
earth's  surface  shown  by  the  hatching.  The  boundary  of  the  central  region 
(within  which  access  is  not  explicitly  forbidden  as  a  result  of  1st  and  2nd  * 

loop  earth  intersections)  la  in  part  formed  by  the  edge  of  the  region 

r 

associated  with  1st  loop  intersections,  and  in  part  by  the  edge  of  the  region 
associated  with  the  2nd  loop  intersections.  As  rigidity  varies  the  region 
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defined  by  the  let  order  shadow  orbits  changes  shape  as  shown  in  Fig.  7. 
The  2nd  order  shadow  lialts  are  wore  stable,  and  change  very  little  over 
this  rigidity  range.  A  change  in  zenith  typically  produces  the  shift  shown 
in  Pig.  8,  and  change  in  asiauth  the  shift  to  the  positions  shown  in  Pig.  9 
and  10. 

The  Band  Shift  Pactor  (Cooke,  1982)  can  be  used  to  quantitatively  express 
the  stability  of  these  structures  (see  Table  1).  It  is  interesting  to  note 
the  relatively  great  stability  of  the  central  second  order  structure,  which 
lies  along  the  geoaagnetlc  equator. 


Table  1.  Band  Stability  Pactor  3R/3P,  where  R  is  rigidity, 
and  P  is  the  the  paraaeter  in  question,  for  the 
points.  A,  B  and  C  in  Plgure  6.  Large  BSP  values 
are  associated  with  relatively  stable  structures. 


LAT 

LON 

ZE 

AZ 

ALT 

GRAZ 

A 

-0.595 

-0.043 

-0.174 

-0.054 

0.015 

-0.018 

B 

-0.750 

-0.030 

0.236 

-0.247 

-0.014 

0.013 

C 

0.109 

-0.004 

-0.245 

0.032 

0.019 

-0.023 

By  examination  of  the  fora  and  position  of  the  shadow  orbit  defined  re¬ 
gions  it  is  possible  to  discount  the  possibility  of  high  zenith  angle  entry 
at  azimuths  well  away  from  the  west,  and  to  rigidities  less  than  a  given  value 
(which  value  depends  on  the  zenith  angle  considered,  e.g.  entry  is  forbidden 
for  rigidity  >  9  GV  at  zenith  angle  of  150*). 


Within  the  central  region  exist  fine  allowed-forbldden  structures  asso¬ 
ciated  with  3rd,  Ath,  and  so  on,  loop  intersections.  The  mapping  of  these 
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structures  has  been  found  Impractical  because  of  the  large  amount  of  computer 
time  required.  In  any  case  the  detailed  form  of  the  structures  depend 
significantly  on  the  choice  of  geomagnetic  field  model  used  in  the  calcula¬ 
tions  (whereas  the  structures  associated  with  1st  and  2nd  loop  intersections 
are  relatively  stable).  It  is  therefore  necessary  to  use  some  other  means 
of  determining  what  access  is  possible  into  this  region. 

Examination  of  the  form  of  the  trajectories  associated  with  approach  to 
points  within  the  central  region  allows  the  question  of  access  from  points 
outside  the  geomagnetic  field  to  be  pursued,  and  permits  general  conclusions 
to  be  drawn  about  whether  access  is  possible  or  impossible  in  given  situa¬ 
tions.  This  examination  has  been  done  by  means  of  sets  of  trajectory  traces 
applying  to  directions  spaced  in  latitude  and  longitude  (typically  at  5° 
intervals  in  latitude,  and  30°  in  longitude).  The  trajectories  are  displayed 
on  a  video  terminal  and  visually  categorized  into  types  (ii),  (ill),  or  (iv). 
Generally  type  (ill)  and  (iv)  orbits  are  readily  distinguished,  although 
there  is  a  less  distinct  difference  between  type  (ill)  and  type  (ii)  orbits. 
The  question  arises  of  when  short  range  quasi-periodicity  turns  into  a  regular 
periodicity.  Not  withstanding  this  point,  meaningful  boundaries  can  be 
established.  Figure  8  shows  the  location  of  the  boundaries  for  this  set  of 
conditions. 

In  general  it  is  found  that  trajectories  originating  close  to  the  geo¬ 
magnetic  equator  show  the  most  simple  and  regular  periodicity  (see  Flge.  3 
and  4).  Further  away  to  either  side  the  trajectories  become  more  complex 
in  their  periodic  form,  and  further  away  again  degenerate  into  unstable  quasl- 
periodlc  or  aperiodic  forms  (see  Fig.  2),  along  which  access  may  occur  from 
outside  the  field  (the  short  u  ,je  shade-  atructures  tend  to  be  responsible 
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for  preventing  access  along  what  otherwise  nay  be  allowed  trajectories).  It 
Is  possible  to  anticipate  whether  trajectories  of  type  (11)  occur  before  the 
shadow  edges  close  off  access,  and  it  Is  therefore  found  unnecessary  to 
actually  locate  an  allowed  trajectory  In  order  to  anticipate  possible  access 
via  trajectories  of  this  type. 

The  offset  of  the  earth's  dipole  gives  rise  to  a  strong  longitudinal 
asymmetry  in  the  form  of  the  regions  (see  Fig.  7  for  example).  There  is  a 
strong  tendency  for  type  (iv)  trajectories  to  be  normally  associated  with 
access  to  points  at  longitudes  close  to  300°,  whereas  type  (ill)  predominate 
over  the  rest  of  the  longitude  range.  Because  the  1st  order  shadow  patterns 
are  mainly  (for  the  boundary  conditions  of  interest)  centered  within  the 
range  90-180°  longitude,  these  regions  can  be  regarded  as  being  populated 
largely  by  type  (iii)  trajectories,  with  type  (ii)  orbits  (if  they  occur 
at  all  in  a  given  situation)  lying  at  points  removed  from  the  geomagnetic 
equator. 

Tn  general,  as  rigidity  decreases  the  size  of  the  mapped  regions  in¬ 
crease.  At  the  same  time,  however,  the  trajectory  configurations  evolve  into 
a  much  more  tightly  bound  periodic  form,  of  types  (ill)  and  (lv).  Particles 
could  conceivably  populate  such  trajectories — possibly  primary  cosmic  rays 
under  time  varying  field  conditions,  or  splash  albedo  particles,  for  example. 

Conclusion 

In  extending  this  preliminary  investigation  it  will  be  necessary  to  make 
an  examination  of  the  zenith,  azimuth,  rigidity,  and  altitude  dependence  of 
of  the  mapped  regions.  In  this  way  a  more  complete  picture  could  be  drawn  up 
of  the  zenith  angle  limits  of  access  -  limits  to  the  regions  containing 
trajectories  associated  with  access  of  the  kinds  possible  and  Impossible, 
probable  and  improbable. 
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Satellite 


LONGITUDE 


Figure  l.  Trajectory  (type  (1))  displaying  short  range  "shadow"  intersection 
with  ataosphere.  The  trajectory,  which  was  continued  on  past  the 
intersection.  Indicated  that  escape  from  field  (or  entry  from 
outside  the  field)  would  otherwise  have  been  possible. 

(Trajectory  arrival  point  parameters  1st.  •  17.5*N.,  long.-  260°  E, 
Ze  -  145*,  As  -  270* ,  Rig.  -  7.5  CV,  Alt.  -  400  km;  Grating  -  30k«; 
geomagnetic  field  IGRF80]. 


L^TuuQg  HEIGHT  ABOVE  GRAZING  ALTITUDE  (Re) 


LATITUDE  HEIGHT  ABOVE  GRAZING  ALTITUDE  (Re) 


iiii] 


LONGITUDE 


Type  (ill)  trajectory  -  showing  regular  quasi-bound  periodic  for* 
intersecting  the  atmosphere. 

{Trajectory  parameters  lat.  7.5,  long.  -  45,  Ze  -  145,  As  -  270 
Rig.  -  7.5,  Alt.  -  400,  Grating  -  30;  geomagnetic  field  IGRF80] 


LATITUDE  HEIGHT  ABOVE  GRAZING  ALTITUDE  (Ro) 
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Figure  5.  Position  of  allowed  and  forbidden  structures  sapped  in  zenith  and 
asiauth  for  the  indicated  rigidity  values,  at  the  point  10s  N 
latitude,  270°  E  longitude  (IGRF80  field).  The  shadow  structure 
limiting  the  high  zenith  angle  entry  in  the  west,  for  the  7.7  GV 
rigidity  value,  is  indicated  by  an  asterlsk(*>. 


Latitude  vs  Longitude 

nith  =  150.000  Azimuth  =  270.000 

titude  =  400.000  Reentry  =  30.000  Year  =  1980.000 


Latitude  vs  Longitude 

>nith  -  146.000  Azimuth  *  270.000  Rigidity  •---  7.500 

ltitude  =  400.000  Reentry  =  30.000  Year  =  1080.000 
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Latitude-longitude  map  for  the  conditions  noted,  showing  categorized  zones  of 
access  (this  Figure  corresponds  to  the  same  set  of  conditions  under  which 
Figure  6  was  produced,  except  that  It  corresponds  to  145°  zenith  rather  than 
the  150°  zenith  of  Figure  6). 


Latitude-longitude  map  drawn  for  same  conditions  as  for  Figure  6,  except  that 
azimuth  is  now  225°,  rather  than  270°  east  of  north. 


Latitude  vs  Loniritude 


Latitude-longitude  map  drawn  for  same  conditions  as  for  Figure  6,  except  that 
azimuth  is  now  315°,  rather  than  270°  east  of  north. 
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